Convergence of the multipole expansion of the polarization interaction 
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The multipole expansion of the polarization interaction between a charged particle and an electri- 
cally charged particle has long been known to be asymptotic in nature, i.e. the multiple expansion 
diverges at any finite distance from the atom. However, it is shown that the multipole expansion of 
the polarization potential of a confined hydrogen atom is absolutely convergent at a distance outside 
the atoms confinement radius. It is likely that the multipole expansion of the dispersion interaction 
of two confined atoms will also be absolutely convergent provided the internuclear separation of the 
two atoms is sufficiency large to exclude any overlap between the electron charge clouds of the two 
atoms. 



O 
Ph. 



> 

o 



oo 

o 



PACS numbers: 34.20. Cf, 31.15.ap 



I. NATURE OF THE POLARIZATION 
EXPANSION 

The long-range interactions between a charged parti- 
cle and an atom is usually described by a polarization 
potential. This interaction is constructed by making a 
multipole expansion of the two-particle coulomb interac- 
tion leading to the following long-range expression for the 
second-order adiabatic polarization potential, 
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In this expression, the ai are the static multipole polar- 
izabilities. The most important term arising from I — 1, 
the static dipole polarizability is sometimes called the po- 
larizability. Equation ([1]) has been long known to be an 
asymptotic expansion. The polarization potential given 
by Eq. ([T]) eventually diverges as I increases at any finite 
value of i? [l|,l3. 

Formal issues of a similar nature also lead to a problem 
in the multipole expansion of the atom-atom dispersion 
interaction. Dalgarno and Lewis 0] originally showed the 
multipole expansion of the H-H second-order dispersion 
interaction. 
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also diverges in the Unsold mean energy approximation 
0, as n — oo at any finite R. Sometime later it was 
shown that the exact multipole series for the H-H inter- 
action was also asymptotic 0, [3 • 

Recently, B-spline basis sets have become very popu- 
lar for solving atomic structure and collisions problems 
Such methods often impose the boundary condi- 
tion that the wave function is zero at some finite radial 
distance, i.e. effectively solving the Schrodinger equation 
in a sphere of finite volume. One aspect of bounding the 



wave function is that the reaction of the wave function 
to an electric field also becomes bounded in space. This 
results in a multipole expansion of the polarization in- 
teraction that can be shown to be absolutely convergent 
provided that the distance from the atom is larger than 
the confining radius of the atom. The imposition of a 
V — oo boundary will result in a multipole expansion 
that is mathematically well behaved at the cost tolerat- 
ing a small error in the size magnitude of the polarization 
potential. 

In this manuscript the Unsold approximation is ap- 
plied to estimate the multipole polarizabilities of a con- 
fined H atom. The resulting polarization expansion is 
then shown to be absolutely convergent when the radius 
is larger than the boundary of the confinement potential. 
This is supplemented by finite radius B-spline calcula- 
tions of the hydrogen polarizabilities for a box size of 
i?o = 14 ao. The results of these explicit calculations 
confirm the Unsold analysis and it is shown that the er- 
ror in the polarization potential arising from the use of 
finite radius polarizabilities does not exceed 3 part in 10^. 



II. POLARIZATION INTERACTION FOR A 
CONFINED HYDROGEN ATOM 

The static multipole polarizabilities of any spherically 
symmetric state are defined 
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where AEgi is the excitation energy from state to state 
i, the sum implicitly includes the continuum, and /g/ is 
the oscillator strength of multipole i connecting the state 
to the excited state i. The f^f are defined [Til. [T^ 
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In this expression C is the spherical tensor of rank £. 
Consider the oscihator strength sum rule 
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This satisfies the identity 
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[13j where N is the number of electrons and ^r"^^') is a 
radial expectation value of the ground state wave func- 
tion. This expression reduces to the well known Thomas- 
Reiche-Kuhn sum rule S'(i)(0) = TV for £ = 1. For the 
hydrogen atom ground state, the sum rule is 
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The exact value of the multipole polarizabilities for the 
hydrogen atom ground state can be written [14i] 
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For a confined atom, the radial expectation value in 
Eq. ^ will be bounded by the radius of the _B-spline 
box, i?0: i-G- 



2£-2 



(9) 



Note that the correctness of Eq. ([6]) for calculations in a 
finite box has been validated to 20 significant figures. 

Eq. ([9]) shows the factorial growth in the sum-rule, 
Eq. ([7]) is eliminated once the wave function is confined. 
For the hydrogen atom ground state one has 
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Confining a hydrogen atom also impacts the energies 
of the states. So the energies used to calculate the oscil- 
lator strengths and the energies used in the energy de- 
nominators of Eq. ([3]). As ^ increases in size, the cen- 
trifugal potential will tend to dominate the coulomb po- 
tential for r < Rq. For £ > \/2ZRo the total poten- 
tial is repulsive and increases monotonically for decreas- 
ing R. When this occurs, the minimum excitation en- 
ergy will be larger than the potential at R = Rq, i.e 
AE = i{e + l)/2Rl - l/i?o + 0.5. For sufficiently large t 
one can write AS « l{£ + l)/2i?g. 

Making the Unsold mean energy approximation Q, 
and setting the mean excitation energy to the minimum 
value, leads to the following large £ approximate expres- 
sion for the hydrogen polarizabilities 
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This ratio is less than one for Ro/R less than 1. Therefore 
the polarization series is absolutely convergent as long as 
R is greater than Rq. 

This absolute convergence has been substantiated by 
calculations of the multipole polarizabilities of the hydro- 
gen ground state. Calculations were performed using a 
B-spline basis with the outer boundary set to Rq = 14 gq 
and results are reported in Table HI The static polariz- 
abilities of hydrogen have been previous calculated using 
B-splines basis sets [l5| but the intent of the present work 
is different from the earlier investigations. 

The radial wave function is written as B-splines 
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where ^'(r) is normalized such that 

|*(r)p dr = 1. 
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The boundary conditions are such that the first B-spline 
is finite at the origin. 
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FIG. 1: The magnitude of the contribution to Vpoi{R ~ 16) 
from each successive multipole of Eq. ([l}. Results for _Ro ~ 14 
ao and the unconfined atom are depicted. 

The polarizabilities listed in Table U show the increas- 
ing difference between the the B-spline calculation of ag 
from the exact calculation at the higher values of £. The 
polarizabilities are written as af" where Rq specifies the 
outer radius of the B-spline box. The B-spline af'° in- 
creases less rapidly than the exact value, a^, and at 

= 20 is 1000 times smaller than the exact polarizabil- 
ity. 

The convergence of the _B-spine polarization potential 
at i? = 16 ao is also shown in Figure[T] The contributions 
of each individual multipole in Eq. ((TJ, i.e. V^^i{R), are 
plotted as a function of i. The divergence in the polariza- 
tion series is apparent as the increments start to increase 
for £ > 15. However, the increments steadily decrease as 
£ increases for the multipole expansion for the B-spline 
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TABLE I: The multipole polarizabilities (in a.u.) of the hydrogen atom ground state. Both exact and B-sphne polarizabihties 
calculated in a sphere of radius -Ro ~ 14 ao are given. The polarization potential at i? = 16 ao is computed with Eq. ([l| 
including terms of successively larger i for both sets of polarizabilities. The column -Eiowost gives the lowest B-spline energies 
for each value of £. 
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a\'^: B-spline 



10-5 X Vri(i?=16) 



10-5 X 1^*1 (i?= 16) 



1 -0.124540597990 

2 -0.043113470410 

3 0.011642665533 

4 0.065428512840 

5 0.122659506743 

6 0.184625743002 

7 0.251827759083 

8 0.324490527950 

9 0.402722627126 

10 0.486577112031 

11 0.576078316293 

12 0.671234886264 

13 0.772046602041 

14 0.878508149199 

15 0.990611293536 

16 1.108346175117 

17 1.231702096001 

18 1.360668006907 

19 1.495232809884 

20 1.635385546034 



4.500000000000000 

15.00000000000000 

131.2500000000000 

2.126250000000000x10^ 

5.457375000000000x10^ 

2.027025000000000x10^ 

1.026181406250000x10** 

6.784199296875000x10^ 

5.671590612187500x10" 

5.846894322018750x10^^ 

7.284255842848359 x 10^^ 

1.078630192114083x10^* 

1.872193833455160x10^" 

3.764357734467175 x 10^^ 

8.679197301530880x10^'' 

2.274460234018827X 10^^ 

6.722293580544535 x 10^^ 



4.499992173038246 

14.99962361522554 

131.2214665586345 

2.123507410263223x10^ 

5.426703060742913x10* 

1.988753009195520x10'' 

9.740932380195944x10'' 

6.020582918170891 xlO** 

4.473394110827833x10^ 

3.838766302059065x10^' 

3.685245972773515 x 10^' 

3.862153915918611x10^' 

4.337864044592644 x 10^' 

5.150165059573718x10^ 

6.397025772360861 xlO^' 

8.248030398257013x10^' 

1.097317979742808x10^' 



2.225432980085533x10^^ 1.499328767223332x10^ 



8.201833248105232x10'^ 



2.096237726358987x10^ 



3.346738528714939 x 10^^ 2.990044876206528 x 10^ 



-3.433227539062500 
-3.477931022644042 
-3.479458973743021 
-3.479555664398503 
-3.479565358643910 
-3.479566765175051 
-3.479567043322079 
-3.479567115152583 
-3.479567138609732 
-3.479567148055899 
-3.479567152652911 
-3.479567155311939 
-3.479567157114797 
-3.479567158530791 
-3.479567159806084 
-3.479567161111560 
-3.479567162618750 
-3.479567164567810 
-3.479567167373771 
-3.479567171846294 



-3.433221567564579 
-3.477923929432083 
-3.479451548358086 
-3.479548114295026 
-3.479557754056124 
-3.479559134030737 
-3.479559398059238 
-3.479559461804639 
-3.479559480306164 
-3.479559486508025 
-3.479559488833742 
-3.479559489785837 
-3.479559490203558 
-3.479559490397286 
-3.479559490491281 
-3.479559490538623 
-3.479559490563226 
-3.479559490576357 
-3.479559490583528 
-3.479559490587524 



V^^^{R = 16) providing a computational demonstration 
of the absolute convergence of the B-spline expansion. 

The imposition of the boundary at Rq ~ 14 ag has re- 
sulted in a polarization potential that is 0.0003% smaller 
in magnitude at i? = 16 oq than the expansion using 
the exact polarizabilities. Most of this difference arises 
from the dipole and quadrupole terms in the multipole 
expansion. 

The demonstrated convergence of the B-spline expan- 
sion has some relevance to the earlier analysis by Brooks 
[iBl- He stated that the origin of the multipole series 
divergence was the use of the expansion 

1 ^ k 

"-K^-g^SrCnnj.cV.) (15) 

beyond its region of validity. This analysis recommended 
that all radial matrix elements involved in the evaluation 
of the polarizability have their radial integrations limited 
to a finite value, e.g. X, which would a lower bound on 
the region for which one writes the polarization potential 
using Eq. ([T]). The interaction, Eq. (fT5|) is effectively 
reduced to 

oo ^ 

^ = E^C'=(^")-C'(X) {r<X) (16) 

where X is a lower bound to the radial coordinate used 
to evaluate V^oi- The resulting multipole expansion was 
then shown to be convergent. Brookes concluded that the 
divergence in the polarization expansion was caused by 
using Eq. (1151) in regions of space where it is not correct. 



This statement is correct, but cannot be used to imply 
that the evaluation of the second-order polarization in- 
teraction using |ri — r2|~^ without making the multipole 
expansion will lead to a convergent expansion in powers 
of . Indeed, it has been shown that by an exact anal- 
ysis in the Unsold approximation that the power series 
expansion R~^ of the hydrogen atom polarization poten- 
tial is asymptotic 1, 2]. However, using a wave function 
that is limited in radial extent, or using an interaction 
that is limited in radial extent 16] will lead to a conver- 
gent expansion of the polarization interaction. 



III. CONCLUSION 

It has been demonstrated that the multipole expansion 
of a confined hydrogen atom polarization potential leads 
to an inverse power series that is absolutely convergent 
provided the distance from the nucleus is larger than the 
box radius. While the result is rigorous for hydrogen, one 
can reasonably assert that the result can be applied to 
any atom. This result has implications that go beyond 
the second-order perturbation expansion of the polariza- 
tion interaction. For example, the Stark effect is known 
to be divergent with respect to order of perturbation the- 
ory [l3]. It is interesting to speculate whether bounding 
the wave function will also remove this divergence. 

The multipole expansion of the atom-atom dispersion 
interaction is also known to be asymptotic in nature 
0, 0, [H| ■ It is reasonable to speculate that this expan- 
sion would also be absolutely convergent for two confined 
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atoms provided the internuclear distance was larger than 
the sum of the two box sizes. It is interesting to note 
that Jansen had suggested that the Vdisp multipole 
expansion would converge provided the basis set used in 
the evaluation of the oscillator str eng th sum rule was 
restricted in dimension and range [la]. The numerical 
results associated with the work of Jansen were rather re- 
stricted in scope. Koide had devised a momentum space 
treatment that resulted in each dispersion parameter be- 
ing multiplied by a radial cutoff of factor [1^ . Although 
it could be shown that the resulting dispersion interac- 
tion was convergent at any finite i? as rt — >■ oo the method 
has only ever been applied to a few terms of the H-H dis- 
persion interaction. 

One aspect of the present analysis should be empha- 
sised. Confining the atomic wave function to a finite 
region of space does constitute the imposition of an arti- 
ficial constraint upon the wave function. However, such 



constraints are routinely applied in calculations of atomic 
structure based on B-spline basis sets. These descriptions 
of atomic structure will automatically have no formal is- 
sues with respect to the convergence of the multipole 
expansion of the polarization interaction. 
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